Slow Dynamics of the High Density Gaussian Core Model 
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We numerically study crystal nucleation and glassy slow dynamics of the one-component Gaussian 
core model (GCM) at high densities. The nucleation rate at a fixed supersaturation is found to 
decrease as the density increases. At very high densities, the nucleation is not observed at all 
in the time window accessed by long molecular dynamics (MD) simulation. Concomitantly, the 
system exhibits typical slow dynamics of the supercooled fluids near the glass transition point. 
We compare the simulation results of the supercooled GCM with the predictions of mode-coupling 
theory (MCT) and find that the agreement between them is better than any other model glassformers 
studied numerically in the past. Furthermore, we find that a violation of the Stokes-Einstein relation 
is weaker and the non-Gaussian parameter is smaller than canonical glassformers. Analysis of the 
probability distribution of the particle displacement clearly reveals that the hopping effect is strongly 
suppressed in the high density GCM. We conclude from these observations that the GCM is more 
amenable to the mean-field picture of the glass transition than other models. This is attributed 
to the long-ranged nature of the interaction potential of the GCM in the high density regime. 
Finally, the intermediate scattering function at small wavevectors is found to decay much faster 
than its self part, indicating that dynamics of the large-scale density fluctuations decouples with 
the shorter-ranged caging motion. 

PACS numbers: 



I. INTRODUCTION 



Essential aspects of the glass transition of the super- 
cooled liquids remain still elusive despite of decades of 
study. Many theories and scenarios have been proposed 
to explain the dramatic slow down of the systems and 
the associated growingcooperative length scales near the 
glass transition point [lMj] . They can explain the exper- 
imental results equally well or equally poorly but none 
of them have been proved to be decisively better than 
other. Even a satisfactory mean-field picture of the glass 
transition has not been established [H, @] ■ Numerical sim- 
ulation of simple model fluids is an ideal route to exam- 
ine the competing theories. Considerable efforts have 
been put forward to gain insight from the dynamical be- 
haviors of simple model glassformers in silico, but com- 
pelling answers are still lacking. There are several reasons 
why the simulation studies are not successful in sorting 
out numerous scenarios and theories. First, the model 
systems are more or less similar; the pair potentials of 
canonical glassformers studied in the past are exclusively 
characterized by short-ranged strong repulsions. Exam- 
ples are Lennard- Jones, its WCA counterpart, soft-core, 
and the hard sphere potentials. Since the strong repul- 
sion dominates thermodynamic and dynamic properties 
of dense fluids, it is hardly surprising that the results 
for these models are qualitatively similar 0, Q . Studies 
of a completely different class of potential systems may 
potentially diversify our views and perspectives on the 
glass transition within the limited accessible time win- 
dows of the simulations. Secondly, the model systems 
are not clean enough. Even the simplest class of model 
glassformers (with a few exceptions [3, [l(| ) are inevitably 
bidisperse or polydisperse in order to avert the nucleation 



to the crystalline phase [7]. This complicates quantita- 
tive assessment of the simulation results. Finally, we still 
lack a realistic model glassformer which conforms to the 
mean-field picture in finite dimensions. Concept of the 
mean-field scenario of the structural glass transition is 
basically borrowed from the mean-field theory developed 
in the spin, glass communities [!, 0, ITTI ? ] . The replica 
theory [l^, [H| and mode-coupling theory (MCT) [lij 
are believed to be the static and dynamic versions of the 
mean-field theory of the glass transition, simply because 
of their apparent resemblance to the spin-glass counter- 
parts. The mosaic pictures of the random first order 
transition theory has been developed as the finite di- 
mension version of this mean field pictures [H, [Til . [To] ]. 
Accumulated simulation data are not inconsistent qual- 
itatively from the prediction of the mean field theories 
but the quantitative agreement between simulation re- 
sults and theoretical predictions are far from compelling. 
The best way to verify the mean-field scenario would be 
to take the mean-field limit by either going to higher 
dimensions or making the system's interactions longer- 
ranged. Recently, simulations for four dimensional sys- 
tems have been performed [l(| EH • Results therein hint 
that the dynamic heterogeneities are suppressed com- 
pared with three dimensional systems and agreement 
with MCT moderately improves [lOj. However, consider- 
ing the current computational abilities, it would be hard 
to simulate the system beyond four dimension, whereas 
the upper critical dimension of the glass transition is ar- 
gued to be eight [l7], EH ■ On the other hand, few studies 
have been done for realistic liquids with long-ranged par- 
ticle interactions [l9l - l2ll |. 

The Gaussian core model (GCM) is a candidate to dis- 
pel all of the above-mentioned concerns and could be an 
ideal and clean bench to compare with various glass theo- 
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ries. The GCM consists of the point particles interacting 
with a Gaussian shaped repulsive potential j22Tj3l| : 



v(r) = e exp[— (rfcr)' 



(1) 



where r is the interparticle separation, e and a are the pa- 
rameters which characterize the energy and length scales, 
respectively. The GCM is one of the simplest models of 
the so-called ultrasoft potential systems which are char- 
acterized by the bounded and long-tailed repulsive po- 
tential [32j] • Recently, we have reported that the one- 
component GCM vitrifies at very high densities [H} . The 
GCM or the ultrasoft particles in general have very dis- 
tinct and exotic proper ties both thermodynamically and 
dynamically [22H3ll l33l. l34j . such as the re-entrant melt- 
ing at high densities, negative thermal expansion coef- 
ficient, and anomalous density dependence of the diffu- 
sion coefficient. There are several studies on the glass 
transition of the ultrasoft particles [3514381 ] and it was 
found that they exhibit rich dynamical behaviors differ- 
ent from conventional model glassformers [35l . |36[ . One 
of the advantages to study the glass transition of the ul- 
trasoft particles is that, due to the mild repulsion tail 
of the potential, the density as well as the temperature 
can be used as a parameter to control the system. Ex- 
ploring the wide range of density-temperature param- 
eter space makes it easier to establish various scaling 
laws, to bridge the gaps between temperature-driven or- 
dinary glasses and density-driven colloidal glasses, and 
to help unifying the concepts of the finite-temperature 
glass transition and the zero-temperature jamming tran- 
sition [37ll38|. However, most studies in the past focused 
on the relatively low density regime, where the generic 
nature of the glass transition is not extremely different 
from that of the conventional model glassformers. The 
systems at low densities, including the GCM, also had 
to be either polydisperse or bidisperse in order to avoid 
crystallization. 

The GCM at very high densities is very differ- 
ent [33|. First of all, the system vitrifies without 
poly(bi)dispersity. The nucleation rate systematically 
decreases as the density increases and the system starts 
exhibiting typical slow dynamics observed in supercooled 
fluids near the glass transition point. Furthermore, the 
dynamics is quantitatively well-described by MCT. Es- 
pecially, the MCT nonergodic transition point extracted 
from the simulation unprecedentedly matches with the 
theoretical prediction. Besides, the violation of the 
Stokes-Einstein relation and the amplitude of the non- 
Gaussian parameter, both of which is the manifestation 
of the heterogeneous fluctuations of dynamics, are sup- 
pressed. We conjecture that these facts can be attributed 
to the long-ranged nature of the interaction potential at 
the high densities where particles are overlapped. These 
results suggest that the high density GCM is not only 
one of the cleanest model glassformers in silico, but also 
the closest to the mean-field model. 

In this paper, we present thorough and complete nu- 
merical analysis of the nucleation and glassy dynamics of 



the high-density and one-component GCM. We not only 
present the exhaustive set of the numerical results but 
also provide with the new evidence which bolsters the va- 
lidity of MCT. Detailed analysis of thermodynamic and 
structural properties of the high density GCM, such as 
the phase diagram and the static structure factors are 
discussed in Ref. 34 1. In the previous study [Hj], we 



have attributed the weak violation of the SE relation and 
smaller non-Gaussian parameter to the suppression of the 
dynamic heterogeneities. We provide stronger and more 
direct evidence that intermittent heterogeneous motion 
is suppressed by monitoring the distribution of the parti- 
cle displacement as a function of time. We also evaluate 
the correlation functions of single and collective density 
fluctuations. Surprisingly we find that dynamics of the 
collective density decouple from the single particle den- 
sity at large length scales, where the former relaxes much 
faster than the latter. This is in stark contrast with the 
ordinary model glassformers for which the slow glassy 
dynamics set in over the whole length scales for both 
collective and single particle densities alike. We com- 
pare these simulation results with MCT predictions and 
find that MCT beautifully captures the decoupling of dy- 
namics at the large length scales. However, we also find a 
subtle but noticeable disagreement of MCT from the sim- 
ulation results at intermediate length scales, where the 
nonergodic parameter (the plateau height of the two step 
relaxation in the density correlators) predicted by MCT 
shows a weak shoulder which tends to grow as the density 
increases. This shoulder is reminiscent of those found for 
the d-dimensional hard sphere glasses at large d evalu- 
ated from MCT [H, Q and may be a signal of breakdown 
of MCT at the mean field limit. 

This paper is organized as follows. In Sec. II, we sum- 
marize the simulation method, theoretical background, 
and the setting of the system. The nucleation dynamics 
from fluid to crystalline phase is discussed in Sec. III. In 
Sec. IV, we present all simulation results on various static 
and dynamical observables. Detailed analysis and care- 
ful comparison of the simulation results with the MCT 
predictions are made. Suppression of the dynamic het- 
erogeneities are also discussed. Finally, Sec. IV concludes 
the paper with a summary. 



II. PRELIMINARIES 
A. Simulation Methods 

We investigate the dynamics of the one-component 
GCM using a molecular dynamics (MD) simulation in 
the NVT ensemble with a Nose thermostat. The sys- 
tem is a cubic cell and a periodic boundary condition 
is imposed. A time-reversible integrator, similar to the 
velocity- Verlet method, is used with a potential cut-off at 
r = 5a [39]. Hereafter, a, e/ks, and oim/e) 1 / 2 are taken 
as the units of the length, temperature, and time, respec- 
tively. The time step is fixed at 0.2, which is sufficiently 
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FIG. 1: State points at which MD simulations were per- 
formed (crosses). Squares with solid line and filled circles 
with dotted line are the solid-fluid phase boundary obtained 
numerically by us 33, 34| and Prestipino et al. [2g], respec- 
tively. The melting and freezing lines are indistinguishable at 
this scale. 



small to conserve the Nose Hamiltonian during the long 
simulation runs. We focus on the four densities, p = 0.5, 
1.0, 1.5, and 2.0, and perform the MD simulations for 
various temperatures in the vicinity of the melting tem- 
perature T m . The state points which we performed simu- 
lation are shown in Fig . [J along with the solid-fluid phase 
boundary line I26l l33l l34j . As discussed in detail in the 
previous study [34]], the melting temperature, T m , at the 
high density regime p > 1 obeys an asymptotic scaling 
logT m oc — p 2 / 3 which was originally conjectured by Still- 
inger 22]. For all densities which we study, thermody- 
namically stable crystalline structure is bec [33, [3J] . We 
run the simulations for the total run time always 50 times 
longer than the structural relaxation time. For example, 
the simulation time was t S i m — 10 7 for the lowest tem- 
perature at p = 2.0. This was confirmed to be sufficiently 
long to neglect aging effect. The first half of the simu- 
lation run was used for the equilibration and we used 
the trajectories of the second half for the analysis of the 
stationary dynamics. For each state point, five indepen- 
dent runs are performed and the results are obtained by 
averaging over those trajectories in order to improve the 
statistics. Configurations obtained from the high temper- 
ature simulation were used as the initial configurations. 
The system size is fixed at N = 3456. The simulations for 
N = 2000 and 9826 confirmed that the finite-size effect 
is negligible. 



transition, MCT is commonly expressed as a set of the 
self-consistent nonlinear equations for correlation func- 
tions. These correlation functions are the intermediate 
scattering function (the correlation of the collective den- 
sity), F(k,t) = (5p(k,0)dp(-k,t))/N, where Sp(k,t) is 
the fc-dependent density fluctuation, and the self inter- 
mediate scattering function or the correlation of the sin- 
gle particle density, F s (k,t) = (p s (k,0)p s (~k,t)), where 

p s (k, t) is the density of a single particle. The time evolu- 
tion of F{k, t) is given by the generalized Langevin equa- 
tion 

n- 2 (k)F(k,t) + F(k,t) + [ ds M(k,t-s)F(k,s) = 0,(2) 



where Sl(fc) = \Jk^Tk 2 /mS(k) is the frequency term. 
S(k) = F(k,t = 0) is the static structure factor. M(k,t) 
is the memory kernel which, according to MCT, is ap- 
proximated as 



M(k,t) = 



pS(k) f dq 



2fc 2 



(27T) 



-V?(q,k- q)F(q, t)F(\k - q\,t).(3) 



Here Vg(g,p) = {k-qc(q)+k-pc(p)} /k is the vertex, where 
c(k) = {1 — l/S(k)}/p is the direct correlation function. 
In Eq. ([3]) , we neglect the short time contribution for the 
memory kernel, which does not affect the slow dynam- 
ics. MCT predicts that F(k,t) undergoes the ergodic- 
nonergodic transition at a finite temperature, T c , below 
which lim^oo F(k, t) — F^k) remains finite. -Foo(/c) is 
referred to as the nonergodic parameter. The nonergodic 
parameter and T c can be evaluated by taking t — > oo of 
Eqs. @ and ([3]), which is expressed as 



l-Fooitysik) 



M OG (fc), 



(4) 



where M 00 (k) is the long time limit of the memory kernel. 
As the temperature approaches to T c from above, MCT 
first predicts that F (k, t) exhibits the two-step relaxation 
behavior characterized by a finite plateau and the slow 
structural relaxation. The height of the plateau is iden- 
tical to F OQ (k) at T — T c . The structural relaxation or 
the alpha relaxation time, r Q , increases and eventually 
diverges at T c . MCT predicts that the increase of r a 
is given by a power law r a ~ \T — T c |~ 7 , where 7 is 
a system-dependent parameter which can be evaluated 
from the MCT equation. 

Likewise, the MCT equation for the self intermedi- 
ate scattering function, F s (k,t), is written in the same 
form as Eq. @, but with the frequency term Q s (k) = 
y ksTk 2 /m instead of f2(fc) and the self memory kernel 



B. Mode Coupling Theory 

In this work, we compare our simulation results for dy- 
namics of the high density GCM in the supercooled state 
with the prediction of MCT. In the context of the glass 



M ^= i^Jwr r^ c(9) } F -^ F ^- ^ 

(5) 

instead of M(k,t) in Eq. ©. The MCT equation for 
F s (k,t) undergoes the nonergodic transition exactly at 
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the same temperature, T c , as for F(k,t), at least for 
most model systems studied in the past (see Ref. [4(| 
for exceptions). By taking the small fc-limit of the 
MCT equation for F s (k,t), we can also construct the 
self-consistent equation for the mean square displace- 
ment (R 2 (t)). MCT predicts that the self-diffusion co- 
efficient D = \im t ^oo (R 2 (t)) /6t follows the power law 
D ~ \T — T c | 7 and vanishes at T c . Note that the power 
law exponent 7 is identical with that for r a . In addi- 
tion to the MCT nonergodic transition and power law of 
the transport coefficients, MCT predicts many important 
dynamical properties such as the dynamic scaling known 
as von Schweidler's law at the plateau regime (the beta 
regime) and the time-temperature superposition at the 
alpha relaxation regime [4lj . 

In order to solve the MCT equations, the static struc- 
ture factor, S(k), is required as an input. We used S(k) 
obtained directly from simulations. For the numerical in- 
tegration of Eq. (|5J) and flSJ , we employed equally spaced 
400 grids with the grid spacing Afc = 0.16. 



III. CRYSTALLIZATION 

Ordinary simple atomic fluids nucleate to form crystals 
quickly as the temperature is lowered below the melting 
point. In this section, we analyze the crystal nucleation 
dynamics of the high density GCM and show that the 
nucleation rate systematically decreases as the density 
increases. In order to monitor the crystallization from 
the homogeneous fluid phase, we use the potential energy 
U and the bond order parameter gg [42[ . The bond order 
parameter is defined by 
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1 N 

»=i 



(6) 



where qi (i) is the Z-th bond order parameter of the i-the 
particle defined by 
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FIG. 2: The time dependence of the bond order parameter 
qe and potential energy U of the representative trajectories 
measured from the time when the system is prepared, (a) 
p = 0.5, T = 2.5 x 10~ 3 , (b) p = 1.5, T = 2.6 x 1(T 5 , (c) 
p = 1.5, T = 2.3 x 1(T 5 , and (d) p = 2.0, T = 2.93 x 10~ 6 . 
The short bold line in each figure indicates the time scale of 
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(7) 



Here qi m (i) is the complex bond parameter of the i-th 
particle given by 



N b (i) 



1 

b .7=1 



(8) 



where B4 is the position of the i-th particle, Nb(i) is 
the number of nearest neighbor particles around the i-th 
particle, and Yi m (r) is the spherical harmonic function of 
the degree I and the order m. 

qe is zero in the fluid phase and q$ w 0.5 for a per- 
fect bec crystal 42] . In Fig. O we show qe and U of the 
five representative trajectories as a function of the lapse 



of time measured from the moment when the system is 
prepared. At a relatively low density p — 0.5 and tem- 
perature just below the melting point T — 2.5 x 10~ 3 
(Fig. [2] (a)), one observes that qe's of all five trajectories 
abruptly increase from zero to a finite value and con- 
comitantly U's decrease. These behaviors are the hall- 
mark of the crystal nucleation. This figure shows that 
the nucleation initiates only after the lapse of time sev- 
eral times longer than the structural relaxation time r a 
which is indicated by the short bold lines in the figures 
(the precise definition and compiled data set of r a are 
given in Sec. IFVl. The degree of supersaturation defined 
by A = 1 — T/T m at this state point is 0.43. Next, we 
look at the higher density p = 1.5. Five runs of qe and U 
at T = 2.6 x 10~ 5 are shown in Fig. (b). Despite of the 
deeper supersaturation (A = 0.55) and much longer sim- 
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FIG. 3: The 54-56 correlation map for the configurations 
obtained at the end of all the five simulation runs at (p, T) = 
(1.5,2.3 x 1CT 5 ) (left panel) and (p,T) = (2.0,2.93 x 1CT 6 ) 
(right panel). Four circles are the characteristic distribution 
for the bcc, hep, fee crystal, and fluid phase. 



ulation runs (over 40 t q ) than Fig. [5] (a), and U do not 
show any sign of nucleation. Decreasing the temperature 
further to T = 2.3 x 1(T 5 where A = 0.6 (Fig. [1(c)), one 
eventually observes the crystallization for the two out of 
five trajectories. Note that it takes decades of the struc- 
tural relaxation time (which itself also increases with the 
degree of supersaturation) before the precipitous nucle- 
ation takes place. At even higher density p — 2.0, all 
five trajectories fail to nucleate even at a very low tem- 
perature T = 2.93 x 10~ 6 with the similar degree of the 
supersaturation, A = 0.6, for the whole simulation runs. 

In order to ensure that the nucleated samples are un- 
ambiguously the bcc crystal and that samples which 
failed to nucleate remain in the homogeneous fluid phase, 
we evaluate new parameters which were recently intro- 
duced by Lechner et al. [43| . They have used the two av- 
eraged bond order parameters q~4(i) and qe{i) and demon- 
strated that the correlation map of them improves ability 
to determine the crystalline structures The aver- 

aged bond order parameter is defined by replacing qi m (i) 
in Eq. (J7]) with the averaged value qi m (i) defined by 



1 N b (i) 

Qim(i) = -~— - V qim(k), 



(9) 



where qi m (k) is given by Eq. (J5]) and the sum runs from k 
over all Nb{i) neighbors of the i-th particles, including the 
i-th particle itself (k — in the sum). In Fig.[3J we placed 
all q~4(i) and <Z6(*) (* = 1, 2, • • • , N) in the correlation map 
for the configurations obtained at the end of simulation 
runs of the two state points (p, T) — (1.5, 2.3 x 10 -5 ) and 
(p, T) = (2.0, 2.93 x 10" 6 ). The four circles represent the 
characteristic areas for the bcc, hep, fee crystals, and 
fluid phase @. The results for (p,T) = (1.5,2.3 x 10~ 5 ) 
show that the two trajectories remain in the fluid phase 
whereas the rest formed the bcc crystal. It is clear that 
no other structures are formed in the course of the sim- 
ulations. Note that the results for the three trajecto- 
ries which nucleated slightly deviate from the bcc region, 
which we presume is due to defects or imperfectness of 
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FIG. 4: The radial distribution function g(r) (left panels) 
and the static structure factors S(k) (right panels), (a) and 
(b) are for p = 1.5 at T = 7.0 x 10" 5 (dashed line) and 
T = 2.4 x 10~ 5 (solid line), (c) and (d) are for p = 2.0 at 
T = 7.0 x 10 -6 (dashed line) and T = 2.93 x 10" 6 (solid line). 
The insets of (b) and (d) are the closeup of S(k) at small k's 
in the semilog plot. Dotted lines in (a) and (c) are the bare 
potential v(r). 



the obtained crystalline structures. On the other hand, 
all the five trajectories for (p,T) = (2.0,2.93 x 10~ 6 ) do 
not show any hint of crystal nucleation and the configura- 
tions remain completely disordered. Hereafter, we focus 
on the densities p — 1.5 and 2.0 because the crystal nucle- 
ation is sufficiently slow that canonical glassy dynamics 
are observed. 



IV. GLASSY DYNAMICS 

A. Structural functions 

Before discussing the slow dynamics in the supercooled 
state, we summarize the fluid structures of the high den- 
sity GCM to demonstrate the difference from those of 
conventional model glassformers. In Fig. [H we plot the 
radial distribution functions g(r) and static structure fac- 
tors S(k) of the GCM for p = 1.5 and 2.0 near and 
below the melting temperatures. Both g(r) and S(k) 
show typical behaviors of dense fluids characterized by 
the prominent peaks near the position and the wavevec- 
tor corresponding to the first coordination shell. Their 
peak heights increase as the temperature decreases. As 
density increases from p = 1.5 to 2.0, the peak position 
of g(r) shifts from r = 0.94 to 0.85 and for S(k) from 
k = 7.8 to 8.4. The noticeable feature of the high den- 
sity GCM is that the tail of the potential v(r) stretches 
beyond the first coordination shell, as demonstrated in 
Figs. 2] (a) and (c). This considerable overlap of parti- 
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cles imparts the character of the long-ranged interaction 
systems to the high density GCM. The long-ranged na- 
ture also appears as the anomalously small S(k) at small 
wavevectors. The insets of Figs. 0] (b) and (d) show that 
S(k 0), or the compressibility, is far smaller than the 
other model fluids at compatible supersaturations, im- 
plying that the density fluctuations at large length scales 
are strongly suppressed. This is the common features 
of the long-range interacting systems. A well-known ex- 
ample is the one component classical plasma pBI |. where 
S(k) vanishes at k — !> 0. More detailed analysis of the 
simulation results for the structural functions and com- 
parisons with the predictions of the liquid state theory 
have been reported in Ref. [34[ . 5(fc)'s obtained here are 
used in the MCT analysis discussed below. 



B. Mean square displacement and self intermediate 
scattering function 

In this subsection, we evaluate various dynamic quan- 
tities and observe their slow dynamics, focusing on the 
trajectories which did not crystallize even when deeply 
supercooled. The mean square displacement (R 2 (t)) = 

A^ _1 EiIi<l-Ri(*) - -Ri(0)| 2 }, the self intermediate corre- 
lation function F s (k,t), and the intermediate correlation 
function F(k,t) are evaluated for the densities p = 1.5 
and 2.0. Fig. [5] shows {R 2 (t)} and F s {k,t) at several 
temperatures well below the melting temperature. These 
figures clearly display the canonical behaviors of the su- 
percooled liquids near the glass transition point. For 
p = 1.5, we could not observe the glassy dynamics below 
T = 2.4 x 1CP 5 because the crystallization intervened. 
At p — 2.0, all trajectories did not crystallize down to 
the lowest temperature which we accessed. In Figs. [5] (a) 
and (c), one observes that, as the temperature is low- 
ered, (R 2 (t)} develops the long plateau regimes followed 
by the usual diffusive behaviors (R 2 (t)) oc t at longer 
times. The appearance of the plateau signals the for- 
mation of a cage of a particle surrounded by its neigh- 
bors and is the hallmark of the supercooled fluid near 
the glass transition point. The value of \J (R 2 (t)) at 
the plateau region is a measure of the sizes of the cages. 
They are about y/{R?(t)} »0.17 for p = 1.5 and 0.14 for 
p = 2.0. These values are slightly smaller than the val- 
ues for conventional model glassformers. For example, 
y (R 2 (t)) s» 0.2 for the Kob- Anderson Lennard- Jones 
mixture (KA model) (46J. 

In Fig. [5] (c) and (d), we plot F s (k — k max ,t) for sev- 
eral temperatures, where fc max is the wavevector where 
S(k) show the maximum peak. -^(fcmax, relaxes ex- 
ponentially at high temperatures. As the temperature 
decreases, a plateau with a finite height appears and it 
stretches over longer times as the temperature decreases 
further, while the plateau height remains almost con- 
stant. This two-step relaxation behavior is another hall- 
mark of the slow dynamic near the glass transition point. 
The terminal relaxation following the plateau is called the 



structural or alpha relaxation. We define the structural 
relaxation time r Q by F s (k jnax ,t = r Q ) = e _1 . In Fig. [HI 
we plot Fjffcjiaj;, t) against the time scaled by r Q . The 
result shows that relaxation curves are collapsed at the 
alpha relaxation regime. This is the universal property of 
the glassy systems known as the time-temperature super- 
position (TTS) [41]. Furthermore, the all curves where 
TTS holds are fitted by a stretched exponential function 
e -(t/r Q ) the exponent (3 ps 0.8. This value is com- 

parable with that for the KA model [46j and for the hard 
sphere mixture [471 ]. 

In Fig. [71 the structural relaxation time r Q and the self 
diffusion constant defined by D = lim t _ J . 00 (i? 2 (i))/6t are 
plotted against the inverse temperature. We plotted D _1 
and adjusted its ordinate so that the data collapses with 
r Q at high temperatures. For both densities, p = 1.5 and 
2.0, r a and D^ 1 drastically increase as the temperature 
is lowered. Both data almost collapse to each other for 
the whole temperatures except for the slight deviation at 
the lowest temperature. As we shall discuss later, this 
is the direct reflection of a weak violation of the Stokes- 
Einstein relation. 

So far, all simulation data show no sign of peculiar- 
ity in the slow dynamics of the high density GCM at 
the qualitative level. They are all similar to conventional 
model glassformers. In order to assess the properties of 
the high density GCM more quantitatively, we compare 
the simulation results with the predictions of MCT. For 
this purpose, we solve the MCT equations Eqs. ©-© by 
numerically integrating the equations in a self-consistent 
manner. As inputs, we used S(k) obtained numerically 
in the previous subsection. First, we compute the MCT 
transition temperature T c by solving Eq. (jU) . The results 

are T c (thcory) = 2.66 x 10~ 5 and 3.17 x 10~ 6 for p = 1.5 
and 2.0, respectively. Here, we denote the transition tem- 
perature as j' c ( thoor y) m order to emphasize that they are 
obtained by solving the MCT equations. The exponent 
7 « 2.7 is also obtained from the MCT solutions. 

MCT predicts that both the self-diffusion coefficient 
and the structural relaxation time follow the power law 
D ,T a oc \T — T c |~ 7 with the same parameters 7 and 
T c . We fitted D~ x and r Q obtained by simulation with 
this MCT power law, using T c as a fitting parameter. We 
denote it as Tc Slm \ By plotting D 1 and t~ 7 against T _1 , 
we found that they both vanish at the same temperature 
and we identified T c (sim) = 2.07 x 10~ 5 and 2.68 x 10~ 6 for 
p = 1.5 and 2.0, respectively (see the insets of Fig.[7J. In 
Fig. El we replotted r Q in Fig. [7J using e = 1 - T/T c (sim) 
instead of 1/T. The results for the KA model [46[ are 
also plotted. These data are scaled by a time unit, to, 
defined by a relaxation time at the short time scale, 
F 8 (k max ,t = to) = 0.95. This figure shows that the re- 
laxation times for both the GCM and KA model ride on 
the MCT power law for the range of temperatures which 
the simulation can access. Collapse of the data of two 
systems on a single power law is a reflection that the val- 
ues of 7's of both systems are close (7 ~ 2.5 for the KA 
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FIG. 5: (R 2 {t)} ((a) and (c)) and F s (fc max ,t) ((b) and (d)). The filled circles are simulation results for p = 1.5 and from left 
to right, T x 10 5 = 7, 4, 3, 2.6 and 2.4 (upper panel), and for p = 2.0 and from left to right, T x 10 6 = 10, 7, 5, 4, 3.4, 3.2, 
3 and 2.93 (lower panel). The dashed line in (c) is the mean square displacement of the KA model at T = 0.475 [4(| shifted 
to fit with the GCM's result at the lowest temperature at long times (see text). The dashed lines in (b) and (d) are the MCT 
solutions obtained using the same reduced temperatures, e, as those for the simulation data. 



model [48} ). This figure also demonstrates that e is a 
good parameter to measure the distance from the onset 
of the glassy slow dynamics for different systems. Here- 
after, we refer to e as the reduced temperature. In Fig. [5] 
(c), we plotted the simulation data of (R 2 (t)) for the KA 
model at T = 0.475 by shifting the time unit in such a 
way that the long time diffusive regime collapses with the 
data for the GCM at T = 2.93 x 1CT 6 and p = 2.0 whose 
reduced temperature is about the same. Almost perfect 
collapse of the results for two distinct systems for the 
whole time window, including the short time ballistic be- 
havior and the entry to the plateau regime, suggests that 
the slow diffusive behavior of the high density GCM is 
qualitatively similar to that of canonical glassformers at 
least above Tc Sim \ where our MD simulation can access. 

However, there are two noticeable differences between 
the high density GCM and conventional model glass- 
formers. First, the MCT transition temperature ob- 
tained from fitting the simulation data, Tc Bim , is un- 
precedentedly close to the theoretical prediction yj* 11001 ^ 
for the GCM. The agreement improves as the density in- 
creases. The deviation of yj slm ) from y c ( thcor y) are on ]y 
32 % for p = 1.5 and 20 % for p = 2.0. It is in stark 
contrast with the KA model for which Tc = 0.435 



and y( theory ) = o.92 with the deviation of more than 



100% [49, 50]. The KA model at T t 



(theory) 



is still a high- 



temperature fluid and F s (k, t) decays exponentially with- 
out a sign of two-step relaxation. Contrarily, the GCM at 
yjthoory) a j reac jy i j eg deep m the region where the plateau 
of F s (k,t) is well developed (see Fig. [5] (d)). Consider- 
able deviation of y c ( sim " ) from yj thcoiy ) f or conventional 
model glassformers is known as one of serious drawbacks 
of MCT. These deviations have been attributed to the ef- 
fect of the activated processes in the ragged energy land- 
scapes, which smears out the clear-cut dynamical tran- 
sition [5lTj54| . Second, the MCT parameters T c and 7 
obtained from fitting simulation data for r Q match very 
well with that obtained from the data of D~ 1 . This is also 
in contrast with the model glassformers such as the KA 
model [Ii|[50| and poly(bi)disperse hard spheres |47Ll55j. 



for which yr im ) ( r the transition density p ( c lm> ) and 7 
obtained from fitting the simulation data vary depending 
on the observables (r Q or D^ 1 ) and also on the compo- 
nents (large or small particles components of the binary 
systems). These variances are partly attributed to the 
presence of strong dynamic heterogeneities which decou- 
ple the diffusion from the structural relaxation time, as 
we shall discuss in the next subsection. 
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FIG. 6: Same as Fig. [5] (d) but plotted against t scaled by 
r a . Filled circles are the simulation results for p = 2.0 and 
T x 10 6 = 5, 4, 3.4, 3.2, 3, 2.93 from right to left. The solid 
line is a fit by a stretched exponential function. 



The direct evidence that MCT works better for the 
GCM than any other model glassformers is the remark- 
able agreement of the simulated F s (k,t) with the MCT 
prediction. In Fig.|S](b) and (d), we plotted the solutions 
of MCT for exactly the same reduced temperatures e as 
the simulation data. Only free parameter is the time unit, 
which is determined solely from the short time dynam- 
ics. Long time behaviors of the MCT solution agree very 
well with the simulation results. MCT also correctly pre- 
dicts the exponent of the stretched exponential relaxation 
j3. The agreement is striking given that for other model 
glassformers, e (and sometimes the wavevectors as well) 
needs to be adjusted at each temperature to obtain a rea- 
sonable fit [49l |56| (an exception is the four-dimensional 
system (lf3|). 



10 £ 



10 4 



10 3 



10 2 



10 e 



10 £ 



^ 10 4 



10 3 




10 7 



10 6 



10 E 



10 4 



10 4 2x10 4 3x10 4 4x10 4 5x10 4 



1/7 




10 £ 



10' 



10 



10 E 



6 i 



10 2 ^ 1 1 1 y 10 4 

10 5 2x10 5 3x10 5 4x10 5 

1/7 

FIG. 7: The temperature dependence of the structural re- 
laxation time (hlled circles) and the inverse of the diffusion 
coefficient (empty squares) for (a) p — 1.5 and (b) p — 2.0. 
Inset: r^T 1 '' 7 and D 1 ^ 1 as a function of inverse temperature, 
where 7 is fixed to 2.7. 



C. Intermediate scattering function 

Next, we look at the intermediate scattering function 
F(k, t). For conventional model glassformers, it is known 
that behavior of F(k, t) is qualitatively the same as that 
of F s [k, t), except for the wiggly fc-dependence of the non- 
ergodic parameter for the former, reflecting the wiggly 
profiles of the static structure factor (see the discussion 
below). Contrarily, for the high density GCM, F(k,t) 
and F s (k,t) differ from each other considerably. Fig. [§] 
shows F(k,t) at two wavevectors. Fig. [5] (a) is the re- 
sult at k = /c max (~ 8.4) which is the peak position of 
S(k). There, the relaxation behavior of F(k,t) is very 
similar to that of F s (k,t), suggesting the relaxations of 
both functions at the interparticle length scales are dic- 
tated by the same relaxation mechanism. Fig. (b) is 



the result at k = 6.4, which corresponds to a slightly 
longer length scale than the interparticle distance. The 
relaxation of F(k, t) is very fast and shows no sign of two 
step relaxation. F(k,t) almost fully relaxed at t ~ 10, 
which is much shorter than the onset time of the caging 
where the plateau of (R 2 (t)) appears (see Fig. [S]). The 
quick decays are followed by the phonon-like oscillations 
and very weak tails persisting over the time scale of the 
structural relaxation time. This tail vanishes at smaller 
fc's. This behavior is in sharp contrast with the KA 
model, where the relaxation time at small wavevectors 
is comparable with that at the interparticle distance and 
the plateau heights remains finite down to very small 
wavevectors [57] . These results indicate that, in the high 
density GCM, the large scale density fluctuations are de- 
coupled from the slow structural relaxation processes at 
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FIG. 8: T a /to as a function of the reduced temperature e for 
the GCM and KA model, to is the short-time relaxation time 
defined by F a (k m ^,t ) = 0.95. 



the shorter length scales. 

In order to see this qualitative difference of F(k,t) 
of the GCM more clearly, we plot the fc-dependence of 
the plateau heights, or the nonergodic parameter, Foo(k) 
and F SiOQ (k) together with the MCT predictions obtained 
from Eq. Q. In Fig. [TUl we show F co (fc)/5'(fc) and 
F SyQO (k) at p = 2.0 (filled circles) and the MCT predic- 
tions at the same density (solid lines). It beyond doubt 
demonstrates that MCT beautifully captures the vanish- 
ing plateau and the decoupling between the self and col- 
lective dynamics at small wavevectors. One observes that 
Foo(k)/ S(k) above fc max remains compatible with that of 
F s ,oo(k), while keeping a wiggly behavior characteristic 
of the collective density fluctuations. Absence of slow 
dynamics at small fc's is a consequence of the anomalous 
structural properties inherent in the high density GCM. 
In the previous subsection, we discussed that the static 
structure factor at the small wavevectors, or the com- 
pressibility, is extremely small compared with those of 
ordinary model glassformers. This makes the amplitude 
of the memory kernel at small fc's negligibly small (see 
Eq. ©). Consequently the large scale fluctuations decou- 
ple from the fluctuations at the length scales of the inter- 
particle distance which trigger the glassy slow dynamics. 
We argue that this decoupling between small and long 
length scales should be commonly observed for the sys- 
tems with small compressibilities which are an universal 
feature of the dense and long ranged interaction systems 
including the Coulomb interaction systems as predicted 
in the framework of MCT [Hj]. 

The nonergodic parameters in Fig. [TU] exhibit another 
subtle but noticeable feature which may have relevance to 
fundamental problems of MCT as the mean field descrip- 
tion of the glass transition. Although MCT reproduces 
the overall behaviors of the nonergodic parameters for 
both Foc^k)/ S(k) and F Si00 (k), its prediction systemat- 
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FIG. 9: The intermediate scattering function at (a) k = 8.4 
and (b) k = 6.2. For both panels, p = 2.0 and the tempera- 
tures are, from left to right, T x 10 6 = 10, 7, 5, 4, 3.4, 3.2, 3 
and 2.93. The inset in (b) shows a closeup of the weak and 
long tails of the main panel. 



ically overestimates the simulation results at the inter- 
mediate wavevectors (in the range of, say, 5 < k < 20). 
As shown in Fig. [10] (b), we find that the simulation 
data for F SyOQ (k) is well fitted by a Gaussian function, 
whereas the MCT nonergodic parameter has a small but 
non-negligible shoulder which the Gaussian function can 
not fit. This shoulder is reminiscent of those observed in 
the MCT solution for hard sphere glasses in large spatial 
dimensions [H, H|. There, we have found that the devia- 
tion from the Gaussian function for F Si00 (k) increases as 
the dimension d increases. This observation has lead us 
to conclude that MCT is not rigorously a bona fide mean 
field theory ||. This glitch of MCT which we found in 
one of the mean field limits, i.e., the high d limit, could 
also show up in another mean field limit, that is, the long- 
ranged interaction limit, which is realized in the high 
density limit of the ultrasoft potential systems such as 



10 



1.0 




o.o -i . . 1 

10 20 30 

k 

FIG. 10: Nonergodic parameters for the collective part 
Foo(k) / S(k) (upper panel) and self part F 3t00 (k) (lower panel) 
of the intermediate scattering functions. Filled circles are the 
simulation data and solid lines are the MCT solutions. The 
dotted line in the lower panel is a fit by a Gaussian function. 



the GCM. This may explain the shoulder of the F St00 (k) 
in Fig. [10] (b). Remember that the anomalously small 
S(k) at small k's is also due to the long-ranged inter- 
action. Interestingly, this small S(k) may explain the 
anomalous shoulder of the MCT solution. By artificially 
enhancing the amplitude of S(k) at small fc's by a minute 
amount and plugging the modified S(k) into the MCT 
equation, we find that the nonergodic parameter F 00 (k) 
at small fc's jumps from zero to finite values. At the 
same time, the shoulder of F Sj0o (k) at the intermediate 
wavevectors disappears and MCT's F Sy00 (k) gets closer 
to the simulation results. This observation implies that 
the long range interaction affects the static properties of 



the large length scales, which eventually amplifies the pu- 
tative non-Gaussian behaviors of the MCT solution. A 
subtle interplay between the long and short length fluc- 
tuations may be quite common for the glass or/and jam- 
ming transition: For example, the the hyper-uniformity 
(vanishing S(k) at small fc) and diverging radial distri- 
bution function at the contact length r = a are known to 
be the two facets of a universal character of the jamming 
transition [59|. 

D. Violation of Stokes-Einstein relation 

For many glassformers, the Stokes-Einstein (SE) rela- 
tion D m T/rj, where r/ is the shear viscosity, is violated 
near the glass transition point and the violation is be- 
lieved to be the manifestation of spatially heterogeneous 
dynamics which grows as the temperature is lowered ffjol ] . 
Indeed, MCT can not capture the SE violation due to its 
mean field character. In this section, we show that the 
SE violation for the high density GCM is suppressed. In 
Fig. HU(a), we plot Dr a for p = 1.5 and 2.0 normalized 
by the values at a high temperature, as a function of e. 
Note that Dr a instead of Dr\ has been plotted, because 
r\ and r a are roughly proportional to each other. In the 
same figure, we have also plotted the results for the large 
and small particles for the KA model [46j]. It is obvious 
that the variations of Dr a for the GCM is much weaker 
than that of the KA model. Similar suppression of the 
SE violation was observed in the four-dimensional hard 
sphere system 10]. 

T a was defined by F s (k,T a ) = e~ 1 at fc = fc max . In 
order to study the length scales which are relevant to the 
SE violation, we generalize the structural relaxation time 
to the fc-dependent form, r(fc), defined by F s (fc,r(fc)) = 
e _1 . Note that r Q = r(fc max ). In the small wavevector 
limit, the self intermediate scattering function behaves as 
F s (k,t) = e- DkH . Therefore, r(fc) ~ 1/Dk 2 as fc ^ 0. In 
the opposite limit, the system should behave as an ideal 
gas, so that F.(k,t) = e - fc B™ 2 * 2 /™ Thus, r(fc) cx 1/fc 
as fc — >• oo [6l|. Fig. [TTJ (b) shows Dk 2 r(k) as a func- 
tion of fc for p = 2.0 and several temperatures. Similar 
analysis for the KA model has been done by Flenner et 
al. [50] ■ At a high temperature T = 7.0 x 10~ 6 where 
the two-step relaxation of F s (k,t) is set off (see Fig. [5] 
(d)), Dk 2 r(k) is nearly constant and almost 1 at small 
wavevectors up to k max . It then decreases as fc increases 
further, followed by a turn over to a mildly increasing 
function. The decrease is a reflection of the vanishing 
of the cages at length scales shorter than the interparti- 
cle distance. The increase at larger fc is a crossover to 
the ideal gas limit where Dfc 2 r(fc) cx fc. The qualita- 
tive behavior remains unchanged at T = 4.0 x 10~ 6 , but 
the drop at fc > fcmax is more pronounced, reflecting the 
stronger cage effect at lower temperatures. At the low- 
est temperature T = 2.93 x 10 -6 which corresponds to 
about e » 0.075, the drop at fc > fc max is more dramatic. 
Furthermore, slight positive bump at 3 < fc < fcmax is 
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FIG. 11: (a) Reduced-temperature dependence of Dr a at 
p — 1.5 (diamonds) and 2.0 (squares). The results for the KA 
model are also plotted (triangles) [4a] . All results are normal- 
ized by those at a high temperature (Dr Q ) ro f. (b) Dk 2 r(k) 
for T x 10 6 = 7.0 (circles), 4.0 (diamonds) and 2.93 (squares) 
at p — 2.0. The arrow indicates fc max , the first peak of S(k). 



FIG. 12: (a) The non-Gaussian parameter a(t) for T x 10 6 = 
10, 7, 5, 4, 3.4, 3.2 3 and 2.93 at p = 2.0. (b) The temperature 
dependence of the maximum value of a(t) at p — 1.5 (dia- 
monds) and 2.0 (squares). The results for the large (up tri- 
angles) and small (down triangles) particles of the KA model 
are also plotted [461 ]. 



observed. This deviation corresponds to a weak SE vio- 
lation observed in Fig. [IT] (a). This behavior is noticeably 
different from that for the KA model for which Dk 2 r(k) 
significantly increases as k increase before dropping near 



E. Non Gaussian dynamics 

Another good measure to monitor the extent of 
the departure from the mean field behavior is the 



non-Gaussianity of the dynamics. At high tempera- 
tures, F s (k,t) or its real space expression, G s (r,t) = 
— Ri(0)\ — r)), also known as the van Hove 
function, becomes almost a Gaussian function. However, 
as the temperature is lowered to the supercooled regime, 
these function substantially deviates from the Gaussian. 
This deviation is also considered to be a manifestation of 
dynamic heterogeneities. To quantify this, it is common 
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to introduce the non-Gaussian parameter denned by 



«(*) 



5(i? 2 (t)) : 



1. 



(10) 



where (R 4 {t)) = N~ x Ei(l^(*) ~ ^(0)| 4 )- In Fig. [12] 
(a), we plot a(t) for /? = 2.0 at several temperatures. It 
shows typical behaviors of the supercooled liquids, char- 
acterized by pronounced peaks at t near or slightly before 
T a whose heights increase as the temperature decreases. 
However, the heights of the peaks are considerably lower 
than other model glassformers at the comparable reduced 
temperatures s. Fig. [15] (b) shows the temperature de- 
pendence of the maximum value of the non-Gaussian pa- 
rameter a max for both p = 1.5 and 2.0. The results for 
the KA model are also plotted (46[. Similarly to the re- 
sult for the SE violation, a max of the GCM is far smaller 
than that of the KA model. Furthermore, one observes 
that a max for P = 2.0 is slightly smaller than that for 
p = 1.5. These results suggest that the dynamic hetero- 
geneities are suppressed for the GCM and the suppression 
is stronger at higher densities. This is another collateral 
support that the high density GCM is more "mean-field- 
like" than other glassformers. 

More direct evidence that the dynamics of the high 
density GCM is closer to a Gaussian process and dynamic 
heterogeneities are weaker can be obtained by monitoring 
the probability distribution of the particle displacement 

; t)is related to the 

m 



r, denoted as P(log 10 r; t). 
van Hove function G s (r, t) by 



P(log 10 r;t) = (lnl0)47rr 3 G s (r,i). 



(11) 



If the dynamics is purely a Gaussian process, G s (r, t) also 
becomes a Gaussian function, 
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(12) 



From Eqs. ([TTj) and (fT2|) . P(log 10 r; t) becomes a function 
of solely r/Vfff; 



P(Iog 10 r;t) = (lnl0)47r 



3r 2 



3/2 



-3r 2 /2(fl 2 (t)> 



(13) 



Thus, the shape of P(log 10 r;i) for a Gaussian process 
should be unchanged as t is varied, but only shifted if 
plotted as a function of log 10 r. The peak height should 
be a constant value of 2.13. In Fig. [T3] we plotted the 
simulated P(log 10 r;t) for p = 2.0 at the two tempera- 
tures; T = 7.0 x 10" 6 (e » 1.2) and T = 2.93 x 10" 6 
(e ~ 0.075). The high temperature result in Fig. [T3] (a) 
shows that P(log 10 r;t) is almost given by Eq. ([TBI ; the 
shape of the function is almost Gaussian and the peak 
height remains very close to 2.13 over the long time. On 
the other hand, the low temperature result in Fig. [T3] 
(b) shows that the peak height of the function becomes 
lower and the width becomes slightly larger at t ~ r a . 



FIG. 13: The probability distribution of the logarithm of the 
particle displacements, (a) The results for T = 7.0 x 1CP 6 
and p = 2.0. From left to right, t = 44, 180, 512, 1024, 2896 
and 5792. (b) The results for T = 2.93 x 10~ 6 and p = 2.0. 
From left to right, t = 500, 32000, 181000, 362000, 1448000, 
and 4096000. 



This non- Gaussian behavior at the beta to alpha relax- 
ation time regime is a common properties of P(log 10 r\ t) 
at a mildly supercooled state. Note that, however, the ex- 
tent of the non-Gaussianity shown in Fig. [T3] (b) is much 
weaker than that of other glassformers such as the KA 
model 50]. P(log 10 r;i) for typical model glassformers 
is known to split into the binodal shape at low temper- 
atures, corresponding to a separation of the constituent 
particles into the mobile and immobile ones. This is one 
of the most salient feature of the dynamic heterogeneities. 
The peak of P(log 10 r; t) in Fig.[T3](b) does not show any 
hint to split into the binodal shape. P(log 10 r; t) of the 
KA model at e = 0.08 (T = 0.47 in the LJ unit), a com- 
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parable reduced temperature as that of Fig. [13] (b), is 
completely separated to the two peaks, corresponding to 
the distribution of mobile and immobile particles. The 
decrease of the peak height of P(log 10 r; t) in Fig. Q2] (b) 
is compatible with that of the KA model at much hi ghe r 
temperature, e = 0.38 (T = 0.6 in the LJ unit) [50]. 
Above results strongly suggest that the dynamics of the 
high density GCM is more Gaussian-like than that of the 
conventional model glassformers and the dynamic het- 
erogeneities are strongly suppressed. 

V. SUMMARY AND OUTLOOK 

In this paper, we presented the detailed analysis of 
dynamics of the high density GCM. The results are sum- 
marized below. 

(i) The crystal nucleation becomes slower as the den- 
sity increases. Analysis of the two orientational bond or- 
der parameters, and q~Q, reveals that the crystal struc- 
ture is bcc at all densities beyond the reentrant point. 

(ii) The system which failed to crystallize shows clear 
two-step and stretched exponential relaxation in the 
(both self and collective) intermediate scattering func- 
tions, which is the hallmarks of the supercooled fluid 
near the glass transition point. All dynamical proper- 
ties which we have analyzed are well described by MCT. 
First, the temperature dependence of the diffusion coef- 
ficient and the structural relaxation time is well fitted 
by the MCT power law. The parameter Tc Sim ^ used to 
fit the simulation data is unprecedentedly close to the 
theoretical prediction. The time dependence of the self 
intermediate scattering function F s (k,t) is well fitted by 
MCT, using the reduced temperature e as a sole pa- 
rameter. Furthermore, the nonergodic parameters for 
both collective and self intermediate scattering functions, 
F oa (k) and F SyOQ (k), are well described MCT. Here we 
find two noticeable differences from the typical glassform- 
ers. First, the shape of F OQ (k) is qualitatively different 
from F S:00 (k) at small wavevectors regime, ~whereF(k,t) 
decays very fast and the nonergodic parameter vanishes, 
whereas F s (k,t) decays very slowly and its nonergodic 
parameter remains finite down to k = 0. We conjec- 
ture that this decoupling of the collective density dy- 
namics from the single particle dynamics is universal for 
the systems with the long-ranged interactions. This indi- 
cates that the large-scale density fluctuation is decoupled 
to the slow structural relaxation processes. Similar de- 
coupling has been predicted from the MCT analysis of 
the systems with the power law interactions v(r) ~ l/r™ 
with small n (58[- Second, the agreement between MCT 
and simulation for F Sj00 (k) is satisfactory but conceiv- 
ably worse than those for other model glassformers such 
as the KA model [13, [E3]- We found a weak shoulder 
at the intermediate wavevectors. This shoulder is remi- 



niscent of those found in the MCT analysis of the hard 
sphere glasses at the high dimensions @. We conjecture 
that the anomalous shoulders are the deficiency of MCT 
which appears only at the mean-field limit. 

(iii) Dynamic heterogeneities are suppressed in the 
high density GCM. The SE violation is very weak and 
the peak height of the non-Gaussian parameter is much 
lower than the conventional model glassformers at the 
comparable reduced temperatures. The weak dynamic 
heterogeneities of the high density GCM was most obvi- 
ous from the observation of the probability distribution of 
the particle displacement P(log 10 r; t). We find no obvi- 
ous change in the shape of P(log 10 r; t) which remains al- 
most Gaussian, though the width slightly widens around 
the beta to alpha relaxation time regimes. Even at the 
lowest reduced temperature, at which the typical model 
glassformers exhibit the very clear binodal distribution 
of mobile and immobile particles, due to the growing dy- 
namic heterogeneities, the probability distribution of the 
GCM remains to be a single peak function. 

We conclude that the high density GCM is the ideal 
model system to study the glass transition. It is not only 
the cleanest glass model in that it is the one-component 
system. But it is also the closest to the "mean-field" 
model in that dynamic heterogeneities are strongly sup- 
pressed and the way how MCT predicts simulation results 
is synchronized with the way it does for the high dimen- 
sional systems. The mean-field nature comes from the 
long-range nature of the interaction potential, which is 
caused by the overlapping of the particles at the high den- 
sities. Both the excellent agreement with MCT and small 
deviation from MCT (the shoulder of F SyOQ (k)) also lead 
us to reconsider the validity of MCT as the the mean field 
theory of the glass transition. Mean-field models of the 
glass transition have been proposed and analyzed by tak- 
ing the long-range limit of the interactions [20| but it has 
never been realized in the simulation box. The another 
mean field limit, i.e., the high dimension limit, is another 
interesting challenge but given the current CPU power, 
going beyond d = 5 would be unrealistic. In this sense, 
the high density GCM might be the first realistic fluid 
model which may be able to bridge the gap between the 
finite dimensional system with the mean-field limit. It is 
tempting to consider the high density limit of the GCM 
where the small parameter 1/p may make the analyti- 
cal treatment of especially the static/thermodynamic pa- 
rameters tractable and leads us the exact mode-coupling 
theory (or alike). 
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